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Tensor network states are used to approximate ground states of local Hamiltonians on a lattice 
in D spatial dimensions. Different types of tensor network states can be seen to generate different 
geometries. Matrix product states (MPS) in D = 1 dimensions, as well as projected entangled 
pair states (PEPS) in D > 1 dimensions, reproduce the D-dimensional physical geometry of the 
lattice model; in contrast, the multi-scale entanglement renormalization ansatz (MERA) generates a 
(D + l)-dimensional holographic geometry. Here we focus on homogeneous tensor networks, where all 
the tensors in the network are copies of the same tensor, and argue that certain structural properties 
of the resulting many-body states are preconditioned by the geometry of the tensor network and 
are therefore largely independent of the choice of variational parameters. Indeed, the asymptotic 
decay of correlations in homogeneous MPS and MERA for D — 1 systems is seen to be determined 
by the structure of geodesies in the physical and holographic geometries, respectively; whereas 
the asymptotic scaling of entanglement entropy is seen to always obey a simple boundary law - 
that is, again in the relevant geometry. This geometrical interpretation offers a simple and unifying 
framework to understand the structural properties of, and helps clarify the relation between, different 
tensor network states. In addition, it has recently motivated the branching MERA, a generalization 
of the MERA capable of reproducing violations of the entropic boundary law in D > 1 dimensions. 

PACS numbers: 03.67.-a, 03.65.Ud, 02.70.-c, 05.30.Fk 



I. INTRODUCTION 

In recent years tensor network states*— have emerged 
as an important theoretical tool to investigate quan- 
tum many-body systems. They offer a novel conceptual 
framework to describe and classify the possible phases of 
matter^l At the same time, as variational ansatze, 
tensor network states are the basis of numerical ap- 
proaches to quantum many-body problems. 

The simplest and best known tensor network state 
is the matrix product state (MPS)i"-l The MPS is at 
the core of the extraordinary success of White's den- 
sity matrix renormalization group (DMRG)£~— , which for 
almost twenty years has dominated numerical research 
in one dimensional lattice models, such as quantum 
spin chains, providing extremely accurate ground state 
properties. The MPS is also used to simulate dynam- 
ics with the time-evolving block decimation (TEBD) 9 -"— 
algorithm and variations thereof, often referred to as 
time-dependent DMRG. Other tensor network states for 
one-dimensional systems include the tree tensor network 
(TTN) 13 ' 14 and the multi-scale entanglement renormal- 
ization ansatz (MERA) 1 ^— , with the later being partic- 
ularly successful at describing ground states at quantum 
critical points. 

Each of the above tensor network states for D = 
1 dimensional systems has a natural generalization in 
D > 1 dimensions. The projected entangled-pair state 
(PEPS)2Ir-2£ generalizes the MPS, whereas D > 1 ver- 
sions of TTN2I>22 and MERA 3 ^ also exist. Among 
those generalizations, PEPS and MERA stand out for 
offering efficient representations of many-body wave func- 
tions, thus leading to scalable simulations in D > 1 



dimensions; and, importantly, for also being able to 
address systems that are beyond the reach of quan- 
tum Monte Carlo approaches due to the so-called sign 
problem, including frustrated spin s 30 ' 39 and interacting 
fermiona 4 ^— . 

An attractive feature of tensor network states is that 
they are largely unbiased variational ansatze, in the sense 
that they are capable of representing many different types 
of ground states through a proper choice of variational 
parameters, as clearly witnessed by two decades of MPS 
explorations with DMRG 5 - - — . By increasing the bond di- 
mension x °f the MPS 4,9 , which governs the size of its 
tensors and therefore the number of variational parame- 
ters, more entanglement can be reproduced and a more 
accurate approximation to the ground state of a lattice 
model is obtained. 

As a matter of fact, an MPS can exactly reproduce any 
many-body state of the system provided that the bond 
dimension \ is large enough 4 ^, although this will typi- 
cally require a prohibitively large value of x, namely a 
value exponentially large in the system size - leading to 
an inefficient representation. What makes the MPS in- 
teresting is that some moderately small value of x is often 
already capable of accurately approximating the ground 
state of a local Hamiltonian in D = 1 dimensions 5 - - — , as 
recently clarified by Hasting s 55 ' 56 . For instance, for most 
gapped Hamiltonians, an accurate MPS approximation 
is obtained already with some finite bond dimension x 
that depends on H but is essentially independent of the 
size of the system. For critical (and thus gapless) systems 
much of the same is true of the MERA-i 5 . - — , which with 
a fixed bond dimension x is capable of accurately repro- 
ducing large scale properties of the ground state, such 
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as the asymptotic scaling of two-point correlators and of 
entanglement entropy. 

Tensor network states in D > 1 dimensions, such as 
PEPS2i^£ and MERA22r— , are also thought to be ca- 
pable of accurately describing a large variety of ground 
states. This is supported both by growing numerical evi- 
dence and by the existence of analytical MER A 35 i 38 and 
constructions for some topologically ordered 
ground states in D — 2 dimensions. However, a sharp 
increase of computational costs with bond dimension x 
implies that, in practice, simulations in D > 1 dimensions 
are restricted to small values of X- This restriction im- 
plies favouring low entangled states over more robustly 
entangled states. And as a result, and unless specific 
preventive measures are taken, tensor network states in 
D > 1 dimensions may artificially favour local order over 
topological order, or may effectively open a gap in a crit- 
ical system. 

It is therefore important to understand how a finite 
bond dimension preconditions the properties of a given 
tensor network ansatz. The goal of our review paper is to 
collect together a number of previous results in this direc- 
tion, which can be found scattered through the literature, 
and to present them under a simple and unifying frame- 
work. We consider homogeneous tensor network states 
with a finite bond dimension x an d look into those struc- 
tural properties (in practice, scaling of correlators and of 
entanglement entropy) that follow from the way the ten- 
sors are connected into a network - and are therefore in- 
dependent of the choice of variational parameters. Those 
properties can be directly associated with properties of 
the (discrete) geometry reproduced by the network. 

A main merit of this presentation is that it exposes, in 
very simple geometric terms, the main structural differ- 
ences between the MPS and the scale invariant MERA 
in D = 1 dimensions, while also explaining why these 
two tensor network states are a natural ansatz to rep- 
resent ground states of gapped and gapless Hamiltoni- 
ans, respectively. The situation turns out to be much 
less clear-cut in D > 1 dimensions, where not only the 
scale invariant MERA but also PEPS can describe certain 
type of gapless ground states; in addition, both PEPS 
and MERA fail to describe another type of more ro- 
bustly entangled, gapless ground states. A second merit 
of our presentation is that it allows us to express, again 
in simple geometric terms, the limitations experienced by 
PEPS and MERA in D > 1 dimensions, while preparing 
the stage to announce how to overcome these limitations, 
namely by using tensor networks corresponding to more 
sophisticated geometries, as recently proposed in Ref. [59|. 

The rest of the paper is divided into sections as fol- 
lows. Section [IT] briefly reviews the typical behaviour of 
correlations and entanglement entropy in ground states 
of local Hamiltonians. In Sect. IIIII we argue that a key 
difference between MPS/PEPS and MERA is in the ge- 
ometry that these tensor network states describe. MPS 
and PEPS reproduce the D-dimensional physical geome- 
try of the lattice model, whereas the MERA describes a 



(D + l)-dimensional holographic geometry, with the ad- 
ditional dimension corresponding to length scale. Then 
Sect. IIVI points out that the decay of correlations in the 
MPS and the scale invariant MERA can be regarded as 
following from the structure of geodesies in the physical 
and holographic geometries, respectively. Similarly, Sect. 
fVl argues that the scaling of entanglement entropy in the 
MPS and the scale invariant MERA can be understood 
to follow from a simple boundary law in some appropriate 
region of the physical and holographic geometries. Sect. 
IVII considers the holographic geometry of ground states 
of gapped systems and a related tensor network state, 
the finite range MERA, which in some sense interpolates 
between the MPS and the scale invariant MERA. 

Finally, Sect. I VIII argues that thinking about tensor 
network states in terms of geometry offers more than 
just an attractive way of presenting previously known 
results. As demonstrated by the recent proposal of the 
branching MERA, it can also stimulate further progress 
in the field. Specifically, the restriction experienced by 
PEPS and MERA to obeying the boundary law for en- 
tanglement entropy in D > 1 dimensions is overcome by 
considering tensor network states that reproduce aD+l 
holographic geometry with a more sophisticated struc- 
ture, including branching in the length scale direction. 



II. GROUND STATES OF LOCAL 
HAMILTONIANS 

Let C denote a lattice in D spatial dimensions made of 
N sites, where each site is described by a complex vec- 
tor space V, and let H : Y® N -> Y® N , with i?t = H, 
be a Hamiltonian that decomposes as a sum of terms, 
with each term involving only a few neighbouring sites. 
We refer to any such Hamiltonian, made of short-range 
interactions, as a local Hamiltonian. The ground state 
I^gs) £ V®^ of H is the state that minimizes the expec- 
tation value (\I>|-ff|^}. Typically, the ground state \^ GS ) 
of a given local Hamiltonian H has a number of struc- 
tural properties that are common to most ground states 
of local Hamiltonians. Here we consider two such struc- 
tural properties: the behaviour of two-point correlation 
functions and the scaling of entanglement entropy. 



A. Correlations 

Let us first consider the two-point correlator 

C( Xl ,x 2 ) = (P X1 Q X2 ) - (P Xl )(Qx 2 ), (1) 

where P Xl and Q X2 denote two local operators acting 
on (a neighborhood of) sites £1,22 € C, and (O) stand 
for the ground state expectation value (vPgsIOI^gs)- In 
most ground states of local Hamiltonians, C{x\,X2) de- 
cays with the distance between positions x\ and X2- In 
the limit of large distances, this happens in one of two 
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characteristic ways. When the Hamiltonian H is gapped, 
correlations decay exponentially^, 

C(x 1 ,x 2 )ne-^- x *Vt, (2) 

where £ > denotes a correlation length. Instead, 
when the Hamiltonian H is gapless, correlations decay 
polynomially^i, 

C(xi,x 2 ) » |xi - x 2 \~ q , (3) 

where g > is some exponent. In this paper we refer to 
gapped systems also as non-critical systems, and to gap- 
less systems as critical systems. A gapped/non-critical 
system that is close to a critical point may display two- 
point correlations that decay according to the more re- 
fined, combined forrn^ 

p -\xt-% 2 \/£ 



which is dominated by the polynomial decay for distances 
smaller than the correlation length, \x\ — x 2 \ <C £, and 
reduces to the exponential decay of Eq. [2] at larger dis- 
tances, \xi — x 2 \ ^> £• 

B. Entanglement entropy 

The amount of entanglement between a region A of the 
lattice and the rest of the system can be measured by the 
entanglement entropy 

S(A) = -tr(p A log 2 p A ), (5) 

where pa is the reduced density matrix for region A, 
obtained from the ground state I^gs) by tracing out the 
rest of the system, denoted B, 

p A =tr B |* GS ){* GS |. (6) 

In D spatial dimensions, most ground states of local 
Hamiltonians obey a boundary law (often also referred 
to as "area law") for entanglement entropy^—, in the 
sense that the entanglement entropy of a hyperblock A 
of L D sites scales as the size \dA\ of its boundary cL4, 

S(A) « \8A\ « L - 1 , (boundary law) (7) 

instead of scaling as the size \A\ = L D of the bulk of 
the block A. However, there are also ground states that 
display logarithmic corrections to the above boundary 

S(A) « L D ~ 1 log 2 (L). (8) 

Specifically, in D = 1 dimensions, gapped systems 
obey a boundary law, which in this case means that the 
entanglement entropy saturates to a constant Sq as a 
function of the block size L 

S(A) < S . (9) 



Instead, critical systems display a logarithmic correction 
to the boundary law, 

5(A)«|log 2 (£), (10) 

where c is the central charge of the corresponding CFT. 
Near criticality, the entropy grows with L as in Eq. [TU] 
until the saturation constant Sq is reached, with the later 
scaling with the correlation length £ as 

So » | log a (0- (11) 

The scaling of Eqs. IMTTlwere first presented in Ref. [l3 
for specific quantum spin chains and are consistent with 
previous entropy calculations in quantum field theory by 
Holzhey et al^i, Callan and WilczekZ2 and Fiola et al^. 
Subsequently, Jin and Korepin 84 formalized the quantum 
spin chain results, whereas Cardy and Calabres o 85 ' 86 for- 
malised and generalized the quantum field theory calcu- 
lations. 

In D > 1, the relation between the scaling of entan- 
glement entropy and the existence of a gap in H is less 
clear-cut. For instance, the study of possible scalings of 
entanglement entropy in the ground state of systems of 
free fermions shows that the boundary law is obeyed by 
gapped systems, as expected, but also for a class of crit- 
ical systems (namely, systems with a Fermi surface of 
dimension L smaller than D— 1), whereas a second class 
of critical systems (with a Fermi surface of dimension 
r = D — 1) display logarithmic multiplicative corrections 
to the boundary law 7 -^— , see Table |TJ Such logarithmic 
corrections are also believed to be present in other gap- 
less systems in D > 1 dimensions, such as Fermi Liquids 
and spin Bose metalsZ^r— . 
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TABLE I. Scaling of entanglement entropy S(A) of a block A 
made of L D sites for the ground state of free fermion models 
on a D-dimensional lattice and with a L-dimensional Fermi 
surface. 



We emphasize that the above simplified characteriza- 
tions of correlations and entanglement entropy touch only 
on those aspects that are needed for subsequent analysis, 
and ignore a number of other important results. For in- 
stance, there are additive corrections to the boundary law 
for entanglement entropy in topological phases, known 
as topological entanglement entrop y 87 ' 88 . Although both 
MERA^ 8 . and PEPSSS can describe topologically or- 
dered systems and in particular account for the topolog- 
ical entanglement entropy, the latter is only an additive 
correction to the scaling of entanglement entropy and as 
such will play no role in our discussions. 
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FIG. 1. (Color online) A local Hamiltonian on a D- 
dimensional lattice defines a discrete D-dimensional geome- 
try. To each ground state we can attach a (D + l)-dimensional 
geometry, where the additional dimension corresponds to 
length scale. 



III. GEOMETRY OF TENSOR NETWORK 
STATES 

A tensor network state expresses the wave-function of 
a lattice system as a collection of tensors (i.e. multi- 
dimensional arrays of complex coefficients) that are con- 
nected according to a network pattern. We refer to the 
abundant literature for further details 1-50 . In this paper 
we are concerned with the use of a tensor network state 
as an efficient, approximate representation of the ground 
state I^gs) of a local Hamiltonian H on a D-dimensional 
lattice. For concreteness, we consider square (or hypercu- 
bic) lattices, although most considerations can be easily 
generalised to other types of lattices. 

There are two geometries that are relevant to the prob- 
lem of representing a ground state, and most of the ex- 
isting tensor network states can be broadly classified ac- 
cording to which of these two geometries their networks 
reproduce, see Fig. Q] 



A. Physical geometry 

First of all, there is the geometry generated by the pat- 
tern of interactions in H, which we refer to as physical 
geometry. In a D-dimensional lattice C, a short-ranged 
Hamiltonian H connects neighbouring sites of C. That 
is, two sites are close to each other in the physical ge- 
ometry if and only if they are also close in the lattice C. 
Therefore the physical geometry is also D-dimensional 
and essentially equivalent to the lattice £ itself. 

By definition, the physical geometry only depends on 
the pattern of interactions in Hamiltonian H and is insen- 
sitive to any further details, such as whether H is gapped 
or gapless. In particular, the physical geometry is also 
largely independent of specific structural properties of its 
ground state I^gs), e.g. decay of correlations or scaling 
of entanglement entropy. 
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FIG. 2. (Color online) (i) Matrix product state (MPS) 
for the ground state of a local Hamiltonian H in a one- 
dimensional lattice. The tensors are connected according 
to a one-dimensional array, in correspondence with the one- 
dimensional physical geometry dictated by the interactions in 
H. (ii) Projected entangled pair state (PEPS) for the ground 
state of a two-dimensional lattice. The tensors are connected 
into a network that reproduces the two-dimensional physical 
geometry. 



An important class of tensor network states consist 
in collections of tensors connected into a network that 
reproduces the physical geometry. For instance, a MPS 
reproduces the one-dimensional physical geometry of a 
spin chain, whereas PEPS reproduces the D-dimensional 
physical geometry in lattice models in D > 1 spatial 
dimensions, see Fig. [2l If we regard the lattice model 
as a discretization of continuous space, then an infinite 
MPS describes a discretized version of the line, and an 
infinite PEPS in D > I dimensions describes a discretized 
version of a D-dimensional hyperplane. 



B. Holographic geometry 

A second geometry is given by the pattern of entan- 
glement in the ground state "Fes) of H. This pattern 
is naturally described by incorporating an additional di- 
mension to the D-dimensional physical geometry. This 
additional dimension is associated to length scale (equiv- 
alently, energy scale), in the spirit of the holographic 
principle^— . Here we refer to the (D + l)-dimcnsional 
geometry generated in this way by the entanglement in 
the ground state |\&gs) as the holographic geometry of 
I'J'gs), and use the scale parameter z, 



= log, A, 



(12) 



where A is a length scale, to label it. The physical geom- 
etry corresponds to setting z = 0. 

In the MERA, tensors are connected so as to reproduce 
the holographic geometry. For instance, the MERA for 
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FIG. 3. (Color online) Multi-scale entanglement renormaliza- 
tion ansatz (MERA) for the ground state of a local Hamilto- 
nian If in a one-dimensional lattice. The tensors form a two- 
dimensional holographic geometry. The horizontal direction 
reproduces the spatial dimension of the lattice model, whereas 
the vertical direction corresponds to the different length scales 
that are relevant to describing the structure of entanglement 
in the ground state of the system. More generally, the MERA 
for a system in D dimensions spans a holographic geometry 
in D + 1 dimensions. 



a one-dimensional system, Fig. [31 spans two dimensions, 
thus describing a discrete, two-dimensional holographic 
geometry, with tensors labeled by a space coordinate x 
and the scale parameter z. More generally, the MERA for 
the ground state of a D-dimensional system reproduces a 
discrete, {D + l)-dimensional holographic geometry. The 
MERA can be regarded as defining a real space renormal- 
ization group transformation, where the scale parameter 
z labels coarse-grained lattices that offer an effective 
description of the system at length scale A = 2 Z . This 
additional dimension allows the MERA to store, using 
different parts of the network, properties of the ground 
state corresponding to different length scales A. 

In contrast with the physical geometry, which is only 
sensitive to the pattern of interactions in the local Hamil- 
tonian H , the holographic geometry depends also on cer- 
tain structural properties of the ground state, such as 
the existence of a finite correlation length £. In order to 
emphasize the differences between the physical and holo- 
graphic geometries, here and in the next two sections 
we consider the holographic geometry of a critical, scale 
invariant ground state, in which the correlation length £ 
diverges and all length scales are equivalent. Correspond- 
ingly, in these sections we will restrict our considerations 
to the scale invariant MERA. Only later, in Sect. IVI1 
we will also address the case of a gapped Hamiltonian, 
where the correlation length £ is finite, and will consider 
a version of the MERA adequate to that situation. 

The connection between MERA and the holographic 
principle was first explored by Swingle 92 . Specifically, the 
scale invariant MERA 15 '— used to describe the ground 
state of a quantum spin chain at criticality can be un- 
derstood as a discrete realization of the AdS/CFT corre- 
spondence. Indeed, the (scale invariant) ground state of 
the critical chain is a discrete version of the vacuum of 
a 1 + I conformal field theory (CFT), whereas the scale 
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FIG. 4. (Color online) As pointed out by Swingle^, the scale 
invariant MERA for the ground state of a quantum spin chain 
can be interpreted as a discrete realization of the AdS/CFT 
correspondence. The ground state of the one-dimensional lat- 
tice model corresponds to a discrete version of the vaccuum of 
a CFTi+i, whereas the MERA spans a two dimensional ge- 
ometry that corresponds to a discrete version of a time slice 
of AdSa+i- The Figure shows a MERA similar to that of Fig. 
|3l but from another perspective, with the scale parameter z 
as a radial coordinate. 



invariant MERA can be regarded as defining a discrete 
version of a 2 + 1 anti de Sitter (AdS) space, see Fig. 
ID [Notice that since we are describing time-independent 
ground states, the time direction is not captured by the 
tensor network and is irrelevant in the present discus- 
sion] . 

C. Homogeneous tensor networks 

Several structural properties of the states that can 
be represented by the above tensor networks are pre- 
determined by the choice of geometry they reproduce. 
The next two sections review the behaviour of correlators 
and entanglement entropy in MPS, PEPS and MERA, 
and relate them to properties of the appropriate geom- 
etry. For simplicity, we mostly consider homogeneous 
tensor networks, in which all the tensors are copies of a 
single tensor (in the scale invariant MERA, two different 
tensors are necessary). 

We call generic a property that is typically observed 
in a homogeneous tensor network where the coefficients 
of the tensor have been chosen randomly. In the follow- 
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FIG. 5. (Color online) Homogeneous tensor network states 
for the ground state in an infinite lattice in D = 1 spacial 
dimensions, (i) A homogeneous MPS is characterized by a 
single tensor that is repeated infinitely many times through- 
out the tensor network, (ii) A homogeneous scale invariant 
MERA is characterized by two tensors, a disentangler and 
an isometry, repeated throughout the tensor network, which 
consists of infinitely many layers. 



ing we consider generic properties of homogeneous MPS, 
PEPS and MERA with a finite bond dimension \- We 
consider states of an infinite lattice £, see Fig. [5J and 
focus mostly on the asymptotic behaviour of such proper- 
ties, namely in the decay of correlations at large distances 
and scaling of entanglement entropy for large blocks of 
sites. 




FIG. 6. (Color online) (i) Two sites xi and xi of the lat- 
tice £ are connected through several paths within the tensor 
network. The geodesic corresponds to the shortest of such 
paths, where the length of a path is measured e.g. by the 
number of tensors in the path. In the example, the shortest 
path between sites xi and X2 contains 6 tensors, and therefore 
the length of the geodesic is D(xi,X2) = 6. Eq. 1131 relates 
the length of geodesies with the assymptotic decay of corre- 
lations in the tensor network (assuming that correlations are 
predominantly carried by the tensors in the geodesic), (ii) 
Region Q.a of the tensor network that contains (the indices 
corresponding to) region A of the lattice C. The boundary 
8Q.A of region J1a consists of the set of indices connecting 
Q.A with the rest of the tensor network. The number n(A) 
of such indices is interpreted as measure of the size |9fM| of 
the boundary 8£Ia- In the example, n(A) = |90a| = 3. An 
upper bound to the entanglement entropy of a region A of the 
lattice C is given in terms of n(A), Eq. 1201 



IV. CORRELATIONS AND GEODESICS 

The asymptotic decay of correlations has long been 
known to be exponential in an MPS 1-3 and polynomial 
in the scale invariant MER A 16 ' 18 ' 19 . In this section we 
point out that such behaviour is dictated by the struc- 
ture of geodesies in the geometry attached to each of 
these tensor network states. For an MPS, the later is a 
rather straightfoward statement; for the MERA, it was 
first noted by Swing 



A. Geodesies within a tensor network 

Given a tensor network state for the state |\&) of a lat- 
tice £, and two sites of C at positions x\ and X2, we can 
define a notion of distance between these two sites within 
the tensor network as follows. First we notice that the 
two sites are connected by paths within the tensor net- 
work, where each path consists of a list of tensors and 
links/indices connecting the tensors. To any such path, 
we then associate a length, as given by the number of ten- 
sors (or links) in the path. Then the distance D(x±,X2) 



between these sites is defined as the length of the shortest 
path connecting them, see Fig. 

Let C(xi,X2) denote a correlation function between 
positions x\ and X2- It turns out that for both the 
MPS and the scale invariant MERA, the decay of correla- 
tions can be expressed in terms of the distance D(x\, X2) 
within the tensor network, 

C(x l7 X2)~e- aD ^ x *\ (13) 

for some positive constant a. This expression assumes 
that the correlations between the two sites are mostly 
carried through the tensors/links in the geodesic path 
connecting them. It originates in the fact that for both 
the MPS (D = 1 dimensions) and the scale invariant 
MERA (in any dimensions), the correlator C(xi,X2) can 
be obtained by evaluating an expression with the (ap- 
proximate) form 

C(x u x2)^v[-(T) D{XI - X 2) -v R , (14) 

that is, a scalar product involving two vectors vl and 
vr and the D(x±, X2)-th power of some transfer matrix 



T. The eigenvalues of matrix T give rise to the possible 
correlation lengths £ in Eq. [5] for the MPS and the possi- 
ble power laws p in Eq. [3] for the scale invariant MERA. 
Instead of reproducing the original derivation of this re- 
sult for MPSi-£ and MERAi&±&i2, here we will focus on 
the geometrical interpretation of Eq. [13] in terms of the 
structure of geodesies. 



B. Correlations in the MPS 

The MPS reproduces the physical geometry of a lat- 
tice C in D = 1 dimensions, and therefore the induced 
physical distance, 

D pbSB (xx,x 2 ) ~ \x\ - x 2 \, (15) 

is simply proportional to the number of lattice sites be- 
tween positions x\ and x 2 , see Fig. [TJi). Replacing the 
physical distance in Eq. [T5] leads to the following asymp- 
totic expression for the correlators of the MPS, 

C MPS (xi,x 2 ) w e -°JW(»i.*0 « e-Ni-^IA (16) 

for some correlation length £ > 0, which indeed repro- 
duces the exponential decay of correlations characteristic 
of gapped systems, see Eq. [2j 

C. Correlations in the scale invariant MERA 

In the scale invariant MERA, two sites at positions Xi 
and x 2 of the lattice C are connected by a geodesic path 
of length 0(log 2 (a;i — x 2 )), see Fig. giving rise to 

the holographic distance 

D bol {xi,x 2 ) w log 2 (|a;i - x 2 \), (17) 

which is consistent with the structure of geodesies in AdS 
space 92 . Replacing this holographic distance in Eq. 13 
leads to the following asymptotic expression for the cor- 
relators of the scale invariant MERA, 

C MERA (x u x 2 )^e- aD ^ x ^ (18) 
„ e -9log 2 (|*i-*2|) _ | JBl _ a , 2 |-g j ( 1Q ) 

for some exponent q > 0, which reproduces the poly- 
nomial decay of correlators characteristic of critical sys- 
tems, Eq. [3] 

D. Correlations in D > 1 dimensions 

In D > 1 dimensions, parts of the same analysis can 
be conducted again for the PEPS and the scale invariant 
MERA. The physical geometry reproduced by the PEPS 
induces a physical distance D phyB (xi, x 2 ) which, as in Eq. 
1151 is proportional to the distance within the lattice £, 
whereas the scale invariant MERA leads to a holographic 
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FIG. 7. (Color online) (i) In an MPS, two spins at positions 
xi and X2 are connected by a path containing |a;i — za| tensors, 
(ii) In a MERA, the same two spins are connected by a path 
that only has 0(log 2 (|a;i — x^D) tensors, in correspondence 
with geodesies in AdS space. 



distance D hoX {xi, x 2 ) analogous to that of Eq. [T7] It is 
also true that, for a generic choice of tensors in the PEPS 
and MERA, we again recover an asymptotic decay of 
correlation functions C(x\,x 2 ) that is exponential and 
polynomial in \x\ —x 2 \, respectively, see Eqs. [161 and IIUI 

However, we point out that for certain (non-generic) 
choices of variational parameters, PEPS can also display 
polynomial decay of correlations, as is the case e.g. of a 
PEPS built from a critical classical partition function 93 . 
Such (non-generic) behaviour is incompatible with the 
assumption implicit in Eq. 1131 namely that correlations 
are mostly carried by the tensors/links included in the 
geodesic path connecting positions x\ and x 2 . In a PEPS 
with polynomially decaying correlation functions, corre- 
lations between two sites are instead obtained from a sum 
of contributions involving the many different paths con- 
necting the two sites within the D-dimensional network, 
and not just from the geodesic paths. Therefore, Eq. [T3] 
does not hold for critical PEPS. 



V. ENTANGLEMENT ENTROPY AND 
BOUNDARY LAWS 

The scaling of the entanglement entropy S(A) of a re- 
gion A is well understood for a MPS^S, a PEPS*! 3 , and a 
scale invariant MERA^. In this section we review these 
results and re-express them as a simple boundary law for 
a related region D,a in the appropriate geometry. 
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(ii) 




FIG. 8. (Color online) Upper bound for the entropy Sl of 
the reduced density matrix pL of a region A of linear size L: 
(i) In an MPS, the tensors describing a block A of L sites, 
that is region f2^ ys , are connected with the rest of the ten- 
sor network by means of two (that is, a constant number of) 
bonds, n(A) = 2. Therefore an MPS can at most reproduce 
a constant entanglement entropy Sl, fa const., which corre- 
sponds to a (physical) boundary law. (ii) In a PEPS for a 
two-dimensional system, the number n(A) of bond indices 
connecting f2^ hys with the rest of the tensor network is pro- 
portional to the size of the boundary of the square region A, 
n(A) ss L. Therefore the entropy scales at most as Sl ~ L, 
which again is a boundary law. 



A. Entropy as the size of a boundary 

An upper bound for the entanglement entropy S(A) 
in a tensor network state is obtained as follows. First 
we consider splitting the tensor network into two parts, 
Qa and Qb, where Qa contains the open indices corre- 
sponding to all the sites in region A and the other part 
Qb contains the open indices corresponding to region B, 
namely the rest of sites in lattice C. Then we count 
the number of bond indices n(A) that connect regions 
Qa and Qb- Since each bond index can contribute at 
most log 2 (x) to the entropy of pa, we obtain the upper 
bounds (see Fig. g^ii)), 



S (A) <n(A) \og 2 ( X ). 



(20) 



Among all such partitions, the one that minimizes n(A) 
is the one that provides the tightest upper bond to S(A). 
From now one, we use Q a and Qb to refer to this optimal 
partition, and n(A) to denote the corresponding minimal 
number of bond indices. The optimal upper bound of Eq. 
1501 is saturated' for MPS, PEPS and MERA, in the sense 
that plenty of numerical evidence shows that for a generic 
choice of coefficients in a homogeneous tensor network, 
the entanglement entropy scales proportional to n(A), 



S(A) « n(A). 



(21) 



In our discrete geometries, we can think of n(A) as a 
measure of the size \8Qa\ of the boundary OQa of the 



minimally connected region Qa, and therefore interpret 
Eq. [5T]as stating that the entropy S(A) is proportional 
to the size of the boundary of region Qa, 



S(A) « \dQ A \- 



(22) 



This expression, equally valid for MPS, PEPS and 
MERA, allows us to always interpret the scaling of en- 
tanglement entropy as a simple boundary law in the ap- 
propriate geometry. The specific scaling of entanglement 
entropy for each tensor network state is then obtained by 
replacing \8Qa\ in Eq. [52] with its explicit dependence 
on the linear size L of region A, as we do next. 



B. Entanglement entropy in MPS and PEPS 

Consider a hypercubic region A made of L D sites. The 
MPS and the PEPS reproduce the physical geometry. 
Therefore the hypercubic region A of the lattice and the 
region f^ lys of the tensor network that contains the open 
indices corresponding to sites in A are essentially equiv- 
alent. In particular, the size of their boundaries is pro- 
portional, |<9f^ hys | w \dA\. In other words, region Q A hys is 
connected with the rest of the tensor network by a num- 
ber of bond indices n(A) proportional to the size \dA\ of 
the boundary dA of the hypercubic region A itself, see 
Fig. E Since \dA\ w L D ~ l , we obtain, 



n{A) pa \8A\ ps L 



D-l 



(23) 



which implies, together with Eq. [5TJ that the entangle- 
ment entropy of the MPS^ and PEPS22 scales with the 
linear size L according to the boundary law of Eq. [7J 
that is 



Smps{A) 

<Speps(^4) 



D-l 



S 



D-l 



(D = l) 
(D > 1). 



(24) 
(25) 



Here we will refer to such scaling of entanglement entropy 
as a physical boundary law or simply boundary law. 



Entanglement entropy in the scale invariant 
MERA 



Given a D-dimensional hypercubic region A of lattice 
C, the minimally connected region Q A ^ in the scale in- 
variant MERA is (D + l)-dimensional, see Fig. [9l with 
the additional dimension labelled by the scale parameter 
z. The number n(A) of bond indices connecting Q A o1 with 
the rest of the tensor network is the sum of T w log 2 L 
different contributions n z (A), 



T-l 



n(A)^J2n z (A), 



(26) 



z=0 



where each contribution n z (A) corresponds to a different 
length scale A = 2 Z , with z G {0, 1, • • • , T — 1}. As ex- 
plained in Ref. l94l the contribution n z (A) is proportional 
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L 



FIG. 9. (Color online) Upper bound for the entropy Sl of 
the reduced density matrix pL of a region A of L contiguous 
sites. In a MERA for a one-dimensional system, the min- 
imally connected region O^" 1 for region A of the lattice is 
connected with the rest of the tensor network by a number 
n(A) that grows logarithmically with the size of region A, 
n(A) ~ log(L). Therefore the entanglement entropy scales 
at most as Sl ~ log(L), which is a logarithmic violation of 
the (physical) boundary law. A more detailed analysis of the 
scaling is found in Fig. 1101 



to the size of the boundary of a region A z obtained from 
region A by means of z coarse-graining steps, where each 
coarse-graining step divides the linear size of the region 
roughly by two, and where the size of the boundary of 
a region is measured by the number of boundary sites 
included in the region. Let us explicitly perform the sum 
in Eq. [26] It is useful to address D = 1 and D > 1 
separately. 

In D = 1 dimensions, given a region A made of L sites, 
each region A z has a boundary dA z made of two sites, 
so that each contribution n z (A) ss \dA z \ = 2 is constant. 
Therefore n(A), which is made of T ks log 2 (L) constant 
contributions (Fig. [TU|) . grows logarithmically with L, 

n(A) «2T«log(Z). (27) 

Then Eqs. [21] and [27] imply that the entanglement en- 
tropy in the scale invariant MERA in D = 1 dimensions 
grows as2i 

S MERA (A) wlog(L) (D = l), (28) 

which reproduces the scaling characteristic of ground 
states of quantum critical systems, see Eq. 1101 

Instead, in D > 1 dimensions, each region A z is a 
hypercubic block of size « L/2 Z (Fig. [TTj) , and therefore 
the size \dA z \ of its boundary dA z scales with L as 

\9A Z \ « f-J . (29) 

Using again that n z (A) is proportional to \dA z \ we find 
that now the contributions n z (A) to n(A) in Eq. [2U 
depend on L. Their sum leads to 

T 

n(A) » L D ~ iy ^2,2~ z « lP~ x . (30) 

z=0 




«( J 4)*log(Z) 

FIG. 10. (Color online) Scaling of entanglement entropy in 
the MERA in D = 1 dimensions, (i) Caricature of region 
A in the lattice C and of the corresponding region f2^ Dl in 
the MERA (see also Fig. [§]). (ii) The total number n(A) 
of bond indices connecting Q,^ 1 with the rest of the ten- 
sor network is the result of log(L) identical contributions, 
each corresponding to a different length scale or value of z. 
Thus, n(A) ~ Iog 2 (L). As a result, the entropy of region 
A in the MERA for D — 1 dimensions scales at most as 
S(A) ~ \og 2 (L), which is a logarithmic violation of the (phys- 
ical) boundary law. 



Then Eqs. and [3D] imply that the entanglement en- 
tropy in the scale invariant MERA in D > 1 dimensions 
grows as 94 

~ L D ~ l (D>1), (31) 

which reproduces the scaling characteristic of ground 
states in most gapped systems and some gapless systems 
in D > 1 dimensions, see Eq. [7] 

Eq. [55] corresponds to a logarithmic violation of the 
(physical) boundary law in D — 1 dimensions, whereas 
Eq. I3H corresponds to the (physical) boundary law. Here 
we reinterpret both Eq. 1281 and Eq. [31] as a holographic 
boundary law, that is, as a boundary law in the region 
1 of the holographic geometry. 

This geometric interpretation of the scaling of entan- 
glement entropy in the scale invariant MERA is inspired 
by (and can be considered a lattice version of) the results 
of Ref. [95l - l97l where the entanglement entropy of a CFT 
is computed using the holographic principle, by noticing 
that it scales as the size of the boundary of a region f^" 1 
with minimal boundary. We emphasize, however, that 
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n(A) «L 



FIG. 11. (Color online) Scaling of entanglement entropy in 
the MERA in D — 2 dimensions. (A similar analysis applies 
to D > 2 dimensions), (i) Caricature of region A in the lattice 
C and of the corresponding region Q 1 ^ 1 in the MERA. (ii) The 
total number n(A) of bond indices connecting fft 01 with the 
rest of the tensor network is the result of log(L) contributions 
n z (A). Contribution n z (A) corresponds to length scale A = 2 Z 
and is proportional to the size \dA z \ ft! L/2 Z of the boundary 
dA z of a coarse-grained region A z . The sum of contributions 
is dominated by the smallest length scale, z = 0, and is thus 
proportional to L. As a result, the entropy of region A in the 
MERA for D — 2 dimensions scales at most as S(A) ~ L, 
which is a (physical) boundary law. 



while Refs. l95l - [97l discuss the scaling of entanglement 
entropy in the actual ground state of a physical theory, 
our present discussion only concerns the scaling of entan- 
glement entropy in a variational ansatz (which we hope to 
be a good approximate representation of ground states). 

One merit of this geometric interpretation is that it 
motivates a strategy to build tensor network states that 
violate the boundary law also in D > 1, as presented in 
Ref. [H and discussed in Sect. IVII1 



VI. HOLOGRAPHIC GEOMETRY IN GAPPED 
SYSTEMS 

The holographic geometry considered so far in Sects. 
IllUfVl corresponds to scale invariant, critical ground 
states, as described by the scale-invariant MERA. This 
particular scenario has been used there to emphasize the 
differences between physical and holographic geometries, 
which are most evident for critical systems. However, all 
ground states, whether corresponding to a critical system 
or a non-critical one, have a holographic geometry. For 
completeness, in this section we consider the holographic 



geometry of the ground states of gapped systems, which 
was first discussed by Swingle^ 2 -. These ground states 
can be represented by a finite range MERAii - a MERA 
with a finite number of layers of tensors, where tensors 
in different layers are in principle allowed to be different. 

Since the finite range MERA is not a homogeneous ten- 
sor network, its structural properties do not only depend 
on the way the tensors are connected into a network - 
different layers of the MERA may contribute differently 
to, say, correlations and entanglement entropy. However, 
even in this case geometrical considerations alone will al- 
ready allow us to reproduce some of the key properties 
that differentiate the ground states of gapped systems 
from critical ones. In addition, near criticality, where the 
correlation length £ is much larger than the lattice spac- 
ing, we will recover aspects of the scaling of correlations 
and entanglement entropy as a function of £. Finally, in 
the opposite limit -namely when the correlation length £ 
is of the order of the lattice spacing- we will see that the 
holographic geometry reduces to the physical geometry 
Correspondingly, for ground states close to this limit, a 
finite range MERA representation becomes equivalent to 
an MPS/PEPS representation. 

A. Correlation length, finite range MERA and 
truncated holographic geometry 

Let us then consider the ground state I^gs) of a gapped 
Hamiltonian H in D dimensions, in which correlations 
decay exponentially with distance according to Eq. [5] (or 
possibly Eq. |4|) and therefore have a characteristic length 
scale, the correlation length £. 

Notice that if we coarse-grain the lattice according to 
a scheme that maps a block of 2 D sites into one site, 
after one coarse-graining step the correlation length £ has 
shrunk by a factor two, £ — > £' = £/2. By applying more 
coarse-graining steps the correlation length will shrink 
further. In particular, after 

**=log 2 (0 (32) 

coarse-graining steps the correlation length will become 
one (in units of separation between lattice sites), and a 
few additional coarse-graining steps, say a fixed number 
Az (independent of £), will render all two-point corre- 
lators negligible (i.e. smaller than some pre-determined, 
small constant). That is, starting with a ground state 
I^gs) with correlation length £, it takes 

z = z ? + Az«O(log 2 (£)) (33) 

steps of coarse-graining to produce a state with negligible 
two-point correlators. 

Here we will assume that after the zo steps of coarse- 
graining (according to an entanglement renormalization 
schem e 15 ' 20 ) the original ground state I^gs) of the sys- 
tem has been transformed into a state that can be well 
approximated by the product state 

|prod} = |0}®|0)<g>---<8>|0), (34) 
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namely a state with no correlations between different lat- 
tice sites. State |prod) describes the ground state at a 
fixed-point of the RG flow that corresponds to a gapped 
phase without topological order. 

In this case, the ground state I^gs) of an infinite sys- 
tem can be represented by a finite range MERA 17 , which 
is made of just a finite number zq of layers of discntanglers 
and isometries. The finite range MERA is not homoge- 
neous, in that the tensors in different layers are allowed to 
be different, reflecting the fact that the properties of the 
ground state are now different at different length scales. 
Here we will consider a particular choice of finite range 
MERA, where the first layers of tensors correspond 
to those in the scale invariant MERA that describes the 
neighbouring critical point, and the remaining Az layers 
are chosen to minimize the ground state energy of the 
gapped Hamiltonian H. This choice is illustrated in Fig. 
\V2\ i) for a finite range MERA in D = 1 dimensions.*^ 

The holographic geometry attached to the ground state 
of a gapped system is still (D + l)-dimensional, but it is 
truncated in the RG direction, with the scale parameter 
z restricted to values in the interval [0,zo]. This trun- 
cation has an immediate effect on the possible decay of 
correlations and scaling of entanglement entropy that the 
tensor network can reproduce. 



B. Correlators 

Let us consider first a correlator between two sites X\ 
and x 2 such that \x\ — x 2 \ is smaller than the correlation 
length £, \x\ — x 2 \ <C £■ In this case the geodesic connect- 
ing the two sites within the truncated holographic geom- 
etry only runs through length scales z smaller than zq 
and therefore is not affected by the existence of the trun- 
cation at z = zq. As a result, the length of the geodesic 
is still logarithmic in \x± — x 2 \ as in a critical system, Eq. 
[171 see Fig. [ToTi). In addition, since the geodesic only 
runs through the homogeneous, scale-invariant region of 
the finite range MERA, z e [0,Z{], the two-point corre- 
lator is expected to decay polynomially as in a critical 
system, Eq. [H 

On the other hand, when \x\ — x 2 \ is larger than the 
correlation length £, \x\ — x 2 \ S> £, the length of the 
geodesic connecting the two sites within the tensor net- 
work grows proportional to \x\ — Xs\, see Fig. I13f ii), 

D hol (xi 1 x 2 ) » \xi - x 2 \, |a;i-a: 2 |>£ (gapped H) 

That is, for sufficiently large \x\ — x 2 \, distances 
D hoX {xi,X2) in the holographic geometry become propor- 
tional to distances D ph y S (a;i, x 2 ) in the physical geometry. 

It would be tempting to say that, in this second regime, 
the structure of geodesies in the finite range MERA, 
Eq. [35l implies that correlations at large distances de- 
cay exponentially. While it is the case that the finite 
range MERA can approximate exponentially decaying 
correlations 110 , these follow from the use of different ten- 
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FIG. 12. (Color online) Finite range MERA for an infinite 
lattice in D — 1 spatial dimensions, (i) = log 2 (£) layers 
of identical disentanglers and isometries obtained from the 
scale invariant MERA for the critical case are followed by 
some fixed number Az of non-homogeneous layers of tensors 
(where the tensors on different layers are allowed to be dif- 
ferent). In this simple example, z^ = 2 and Az = 1, so that 
the total number of layers of tensors is zq = z^ + Az — 3. 
The finite range MERA represents a ground state |\Pgs) of 
a gapped system that can be transformed into the unentan- 
gled state |prod) after three layers of coarse-graining, (ii) 
The finite range MERA can be combined with another ten- 
sor network (e.g. MPS in the figure) in order to represent 
a ground state |$gs) of a gapped system that flows towards 
an entangled fixed point ground state l^t.p.) under coarse- 
graining transformations. The MPS at the top of the MERA 
can also be used to accurately describe the exponential decay 
of correlations that dominate the limit \xi — x%\ 2> £ of Eq. 
II 



sors in top Az layers of the network, and cannot be inter- 
preted in simple geometric terms. Nevertheless, what is 
clear from geometric arguments (together with some un- 
derlying transfer mechanism for the propagation of corre- 
lations) is that a truncated holographic geometry can no 
longer give rise to polynomially decaying correlations at 
long distances, \x\ — x 2 \ ^> £, since the length of geodesies 
is no longer logarithmic, Eq. 1351 

In conclusion, we have argued that the truncated holo- 
graphic geometry of the ground state of a gapped system 
gives rise to a modified structure of geodesies that is com- 
patible with the decay of two-point correlators expressed 
in Eq. 31 with polynomial decay for |aii — x 2 \ < ^ and 
exponential decay for \xi — x 2 \ 3> £. 
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FIG. 13. (Color online) Geodesies in the finite range MERA. 
(i) When \xi — Xa| is smaller than the correlation length £, 
the geodesic connecting sites xi and X2 within the tensor net- 
work is identical to the scale invariant case, and its length 
is therefore logarithmic in \x\ — x%\. (ii) When \x\ — xi\ is 
larger than the correlation length £, the geodesic connecting 
the two sites sees the presence of the truncation and grows 
proportional to \x\ — xi\. The structure of geodesies in the 
truncated holographic geometry is therefore compatible with 
the refined decay of correlations of Eq. [4] 




L 



FIG. 14. (Color online) Minimally connected regions in the 
finite range MERA. (i) When region A is smaller than the 
correlation length £, L <C £, the minimally connected region 
f}J\ o1 within the tensor network is identical to the scale in- 
variant case, and the size n(A) = | S^^ 1 1 of its boundary Q^ 1 
is therefore logarithmic in L. (ii) When region A is larger 
than the correlation length (, L> the size of the bound- 
ary O^ 1 saturates to a constant (as a function of L) that is 
proportional to the number zq of layers in the finite range 
MERA and thus grows logarithmically with f. The structure 
of minimally connected regions in the truncated holography 
is therefore compatible with a saturation of the entropy, Eq. 
[9] with a saturation value So that scales as log 2 (£), Eq. 1111 



C. Entanglement entropy 

In D — 1 dimensions, the scaling of entanglement en- 
tropy in the finite range MERA (for gapped systems) is 
also different than in the scale invariant MERA (for criti- 
cal systems) , with the difference having a straightforward 
geometric interpretation. 

Let us consider Fig. [TU If we first consider a region A 
with length L smaller than the correlation length £, L -C 
£, then the minimally connected region fi^ 1 in the tensor 
network only involves length scales z smaller than z^ . As 
illustrated in Fig. H~4T i). in this case the size n(A) = 
<9f^ 0l | of the boundary dfl^ 1 scales as log( L) as in the 
scale invariant MERA, Eq. [27J and the entanglement 
entropy grows as in a critical system, Eqs. 1101 

However, when the size L of region A is larger than 
the correlation length £, L 3> £, the minimally connected 
region f2^ 01 in the finite range MERA has a boundary 
dfl^ 1 that saturates to a constant size Idfl^'l, see Fig. 
HKii), with 

\dn h /\ = n(A) « zz = log 2 (0 (D = 1, gapped H) 

(36) 

That is, region f^" 1 is now connected with the rest of the 
tensor network through a number n(A) of bond indices 
proportional to z^ = log 2 (£)> which is independent of L. 
This implies a constant upper bound for the entangle- 
ment entropy 

S MERA (A) « L D - X ^S {D= 1, gapped H) (37) 

and therefore the finite range MERA obeys the bound- 
ary law of Eq. |H1 In addition, the entropy saturates 
to a constant So ~ z^ that grows logarithmically with 
the correlation length thus also producing a scaling 
compatible with Eq. [TT] 

Therefore we see that the structure of minimally con- 
nected regions in the truncated holographic geometry 
reproduces well the scaling of entanglement entropy in 
gapped systems, both for block lengths L larger and 
smaller than the correlation length £. 

In D > 1 dimensions, the truncation of the holographic 
geometry to z < zq due to the presence of a finite corre- 
lation length £ does not alter the scaling of entanglement 
entropy, which is dominated by the z — contribution 
(see Fig. [T5|) . and therefore the finite range MERA still 
obeys a boundary law, 

S MERA (A) a \dn h /\ » L D - 1 (D > 1, gapped H) (38) 



D. Equivalence between holographic and physical 
geometries 

We have seen that for gapped systems, the holographic 
geometry is truncated at a value zq of the scale param- 
eter z corresponding (up to a constant) to = log 2 (£), 
where £ is the correlation length. We have also seen that 
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FIG. 15. (Color online) Scaling of entanglement entropy in 
the finite range MERA in D = 1 and D — 2 dimensions, 
(i) For D = 1, region fix 1 is a truncated version of that in 
Fig. 1101 (ii) Only length scales smaller than £ contribute to 
the total number n(A) of bond indices connecting f^°' with 
the rest of the tensor network. This number is therefore up- 
perbounded by a constant, which grows with the correlation 
length as log 2 (a;i). (iii) For D = 2, the region fi^ is a trun- 
cated version of that in Fig. 1111 (iv) Again, only length scales 
smaller than £ contribute to the total number n(A) of bond 
indices connecting Q,^ 1 with the rest of the tensor network. 
However, this does not change the linear dependence of n(A) 
in L. 
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FIG. 16. (Color online) The finite correlation MERA can 
be converted into an MPS with a sufficiently large, but finite 
bond dimension Xmps as given by Eq. 1391 Specifically, each 
bond index of the MPS has to account for O(log 2 (£)) bond 
indices of the MERA. When £ is small, the two-dimensional 
holographic geometry of I^gs) an d the one-dimensional phys- 
ical geometry of H are essentially equivalent, and the ground 
state can be accurately described by either an MPS or a 
MERA. Figs, (i)-(iii) illustrate in diagrammatical notation 
how to convert a finite range MERA into a MPS, where the 
bond dimension of the MPS grows exponentially with the 
number of layers in the MERA. 



the presence of the truncation in the holographic geome- 



try implies that the length D hol (xx, £2) of geodesies and 
the size [cT^'l of the boundary of minimally connected 
regions in the holographic geometry scale asymptotically 
as in the physical geometry. As a matter of fact, when the 
correlation length £ is of the order of the lattice spacing, 
so that zq is a small number, it is no longer possible to 
distinguish between holographic and physical geometries 
at all. 

Correspondingly, as we discuss below, in D = 1 dimen- 
sions the finite range MERA can be efficiently mapped 
into an MPS. This mapping if still possible for a large 
correlation length £, but the resulting MPS has a bond 
dimension x that grows with £ and diverges at a critical 
point. 



E. From MERA to MPS 

As illustrated in Fig. [16] for gapped systems in D = 
1 dimensions, a finite range MERA made of zq layers 
of tensors and with bond dimension Xmera can be re- 
expressed as an MPS with bond dimension Xmps given 
by^ 



Xmps ~ (Xmera) 



zq~ z ( = log 2 (£). 



(39) 



Therefore, when the correlation length £ is small, the 
holographic geometry is a narrow strip and the MERA 
can be re-expressed as an MPS with a small bond dimen- 
sion Xmps- 

However, as one gets closer to a quantum critical point 
and the correlation length £ becomes larger, the holo- 
graphic geometry becomes a strip with larger width zq. 
The MERA can still be re-expressed as an MPS, but with 
a bond dimension Xmps that grows exponentially with 
the number zo of layers in the MERA, Eq. Since 
the computational cost of MPS algorithms scales as a 
power of Xmps j that is exponentially with zq, numerical 
simulations must be restricted to small values of zq. 

Finally, at the critical point, where £ diverges, the 
holographic geometry extends indefinitely in the coarse- 
graining direction z, and a MERA with finite bond di- 
mension Xmera can no longer be replaced by an MPS 
with finite bond dimension Xmps- 

In D > 1 dimensions, a finite range MERA with a 
finite bond dimension Xmera can also be re-expressed as 
a PEPS of finite bond dimension Xpeps- However, the 
bond dimension Xpeps does not grow significantly with zq 
and, as a matter of fact, even the scale invariant MERA 
in D > 1 (with zq = 00) can be exactly represented by 
a PEPS with finite bond dim ensio n Xpepsj as recently 
shown in Ref. [99| (see also Ref. 1 1 12T ) . 



VII. DISCUSSION 

In this manuscript we have reviewed a number of re- 
sults concerning correlations and entanglement in ten- 
sor network states and presented them in a unified way 
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by pointing out that they can be interpreted as ge- 
ometric properties of some underlying discrete geom- 
etry. Specifically, MPS and PEPS have been argued 
to describe a Z?-dimensional physical geometry dictated 
by the interactions of a local Hamiltonian H in D- 
dimensions, whereas the MERA has been seen to de- 
scribe a (D + l)-dimensional holographic geometry as- 
sociated to the ground state I^gs) of H. Our presen- 
tation clearly emphasizes the main structural differences 
between MPS/PEPS and MERA, and interprets the de- 
cay of correlations and the scaling of entanglement en- 
tropy in terms of geometric concepts such as geodesies 
and regions of minimal surface within the relevant geom- 
etry, as summarized in Eqs. [13] and [22] The geometrical 
interpretation is also the natural language to connect the 
MERA with the holographic principle. 

We conclude the present review with two brief discus- 
sions of related issues. The first is a practical warning 
for future tensor network practitioners. The second is a 
pointer to on-going developments that have been moti- 
vated by the geometric perspective described here. 



A. MPS for critical systems in D — 1 dimensions; 
and for systems in D — 2 dimensions. 

First, a word of caution on the use of geometric consid- 
erations to characterize tensor network states is in order. 
Our discussion has mostly focussed on the asymptotic 
decay of correlations and asymptotic scaling of entan- 
glement entropy. In D = 1 dimensions, this analysis 
pointed at the MPS as a natural representation for the 
ground state of gapped systems, and at the scale invari- 
ant MERA as a natural representation for the ground 
state of critical systems. However, this should not be 
understood as implying that an MPS cannot be used to 
study ground states of critical systems in D = 1 dimen- 
sions, or even ground states in D = 2 dimensional lat- 
tices. 

In a homogeneous MPS with finite bond dimension 
Xmps, two-point correlators are indeed constrained to 
asymptotically decay exponentially, Eq. II 61 whereas 
entanglement entropy must saturate, Eq. [M] thus re- 
producing the scaling of Eqs. [5] and M characteristic of 
gapped systems. However, for some intermediate values 
of distance \xi — x?\ and size L, an MPS can still accu- 
rately approximate a polynomial decay \x\ — X2\~ q of cor- 
relations and a logarithmic growth log 2 (i) of entangle- 
ment entropy. More specifically, a finite bond dimension 
Xmps m the MPS has been seeni22^i22 to introduce an 
artificial, finite correlation length £ x (where £ x depends 
on central charge c of the CFT that describes the critical 
point under consideration) such that the correct scaling 
of correlators and entropy is reproduced for \x\ — X2 and 
L smaller than £ x . 

A relatively mild scaling of computational costs with 
the bond dimension of the MPS, namely as Xmps to the 
third power, implies that very large bond dimensions (of 



the order of thousands) can be afforded with reasonably 
modest computational resources, leading to large values 
of the effective correlation length £ x . This, together with 
the use of finite size scaling techniques, make the MPS a 
very suitable tool to study critical ground states, which 
explains the success of DMRG also for critical systems^—. 

Similarly, an MPS may appear as an unlikely candi- 
date to represent ground states of D = 2 lattice models, 
since the only way it can afford reproducing the bound- 
ary law of entanglement entropy in D = 2 dimensions, 
Eq. [7J is through a bond dimension Xmps that grows ex- 
ponentially in the linear size of the lattice. Once more, 
however, using an MPS with very large Xmps (which can 
again be afforded due to the relatively mild scaling of 
computational costs with Xmps) and finite size scaling 
arguments, an MPS has been successfully used to study 
ground states of two-dimensional lattice models^r-iSi. 



B. Beyond the entropic boundary law in D > 1 
dimensions. 

The present analysis has also reminded us of an impor- 
tant limitation of PEPS and MERA in D > 1 dimensions. 
Recall that these tensor network states are constrained 
to obey a strict boundary law for entanglement entropy, 
Eq. [JJ However, there is an important class of gapless 
systems in D > 1 dimensions whose ground states dis- 
play a logarithmic violation of the boundary law, Eq. 
[U These systems include Fermi gases and liquids with a 
(D — l)-dimensional Fermi surface, as well as spin Bosc- 
metals with an analogous Bose surface^— . How may 
we go about using a tensor network state to represent 
such ground states? 

In the case of PEPS, the boundary law cannot be easily 
overcome, since it is an intrinsic property of the physi- 
cal geometry that the ansatz reproduces. Mimicking the 
previous discussion on the use of MPS to study critical 
systems, one could, perhaps, study ground states with 
a logarithmic violation of the boundary law with by a 
PEPS by considering finite systems and by suitably in- 
creasing the bond dimension Xpeps with the system size. 
Then finite size scaling techniques could be used to ex- 
trapolate finite size results to the thermodynamic limit. 
However, the cost of PEPS simulations grows as a much 
larger power of the bond dimension Xpeps than in the 
case of MPS, confining Xpeps to small values and seri- 
ously limiting the viability of this strategy. 

Similar considerations apply to the MERA: a system- 
atic increase of bond dimension xmera as larger systems 
in D > 1 dimensions are considered seems unviable, due 
to the sharp increase of computational costs with Xmera ■ 
However, remember that in D = 1 dimensions the loga- 
rithmic violation of the boundary law, Eq. 1101 could be 
interpreted as a boundary law in the holographic geom- 
etry, Eqs. [HI and [55J This strongly suggests a possible 
alternative. Indeed, it is natural to wonder whether, once 
more, the logarithmic violations of the boundary law in 
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FIG. 17. (Color online) Holographic branching in a D — 
1 dimensional system, (i) Graphical representation of the 
holographic geometry in the absence of branching. Space is 
labelled by coordinate x and whereas the scale parameter z 
labels the different length scales A = 2 Z in the system, (ii) A 
region A in the D = 1 system defines a minimally connected 
region T^ 01 in the the holographic geometry, as discussed in 
Sect. |V] (iii) Boundary dT^ 1 of the minimally connected 
region T 1 ^ 01 . (iv) Graphical representation of the holographic 
geometry in the presence of branching at the value z — z* of 
the scale parameter, corresponding to length scale A* = 2 Z . 
(v) The same region A of (ii) gives now rise to a minimally 
connected region F^ 01 which is also affected by the branching 
(provided that the size L of region A is larger than the length 
scale A* at which branching occurs), (vi) As a result of the 
branching, the boundary dT\ oX of the minimally connected 
region F^ 01 is larger. A larger boundary leads, together with 
Eq. 1221 to a larger amount of entanglement entropy. 



D > 1 dimensions, Eq. [8l can follow from a boundary law 
in a more elaborated, yet unknown, {D + l)-dimensional 
holographic geometry. A generalized MERA that would 
reproduce this holographic geometry would then auto- 
matically display a logarithmic violation of the boundary 
law. 



C. Holographic branching 

As recently discussed in Ref. [H, it turns out that, 
indeed, one can engineer holographic geometries such 
that Eq. [8j as well as many other forms of scaling, can 
be understood to follow from a holographic boundary 
law. A key ingredient in these holographic geometries is 
the presence of branching, by means of which a single 
(D + l)-dimensional geometry associated to small length 
scales (high energies) becomes two independent (D + 1)- 
dimensional geometries at large length scales (lower en- 
ergies), see Fig. [T7](i). Physically, holographic branching 



(i) (ii) a (iii) 




S(A) « L S(A) « L S(A)*L log(L) 

FIG. 18. (Color online) Schematic representation of differ- 
ent holographic geometries in D — 2 dimensiona l sy stems in 
terms of the holographic tree introduced in Ref. 153 . (i) The 
holographic geometry of a gapped system has a single branch 
with a finite extension in the z direction, namely z £ [0, zo] 
where zo = 0(log(£)) as discussed in Sect. I VII This geometry 
corresponds to a finite range MERA. (ii) The holographic ge- 
ometry of a gapless system that obeys the entropic boundary 
law may also consist of a single branch, but this extends in- 
definitely in the z direction, as in the scale-invariant MERA. 
(iii) Holographic geometry with an infinite number of branch- 
ing points, capable of reproducing the logarithmic violation of 
the boundary law characteristic of the ground state of several 
gapless systems, including free fermions, see Tabled] Notice 
that the holographic geometry of (i) and (ii) allow us to dis- 
tinguished between two types of ground states, corresponding 
to gapped and gapless systems, that obey the boundary law. 
In this sense, the holographic geometry can be used to issue a 
more refined classification of ground states according to their 
pattern of entanglement. 



describes the decoupling of a single theory into two the- 
ories (or sets of degrees of freedom) that do not interact 
with each other at energy scales lower than some decou- 
pling energy - equivalently, at length scales larger than 
some decoupling length A. Thus, at length scales smaller 
than A there is a single lattice model, whereas at length 
scales larger than A, the lattice model breaks into two 
independent lattice models. 

Fig. [T7] illustrates how the presence of holographic 
branching affects the amount of entanglement entropy in 
the ground state. The holographic region associated 
with a physical region A of linear size L larger than A also 
branches into two pieces. As a result, the entropy S(A) 
receives contributions from two pieces of the boundary 
80, a- As discussed in Ref. 591, it turns out that a se- 
quence of holographic branchings occurring at different 
length scales, as represented by a branching tree (see Fig. 
ITSn ii) for an example) leads to a wide range of forms of 
scaling for the entanglement entropy S(A) of a region A 
in the original lattice, including Eq. [SJ The resulting 
tensor network state, the branching MERA, reproduces 
these model elaborated holographic geometries and has 
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been shown to efficiently represent e.g. the ground state 
of a D = 2 dimensional fermionic lattice model with a 
one-dimensional Fermi surface. 

The study of holographic geometries with (possibly 
multiple) branching points opens up a number of excit- 
ing new possibilities, presently under consideration. On 
the one hand, it motivates a revision of the RG flow and 
its structure of fixed points. As we progress towards low 
energies, a single theory may branch into (perhaps in- 
finitely many) other theories. In particular, new fixed 
points of this revised RG flow, including branching at all 
length scales, seem to include certain D = 2 dimensional 
systems with a one-dimensional Fermi surface. 

On the other hand, the holographic geometry also of- 
fers a new venue to characterize entanglement of ground 



states. The pattern of branching (as given by a holo- 
graphic tree) of a ground state I^gs), as well as the ex- 
tent of each branch in the scale direction z, leads to a 
new classification of ground states that subsumes the one 
provided by considering only the scaling of the entangle- 
ment entropy S(A) of a region A of the lattice, see Fig. 
[TBI For instance, while gapped systems and some gapless 
systems in D > 1 dimensions cannot be distinguished by 
the scaling of entanglement entropy (since they all obey 
the boundary law of Eq. [7} , their holographic geometry 
is clearly distinct. 
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rated when all the eigenvalues of pA are equal to l/r(pA)- 
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with another tensor network, to describe other gapped 
phases of matter, such as (possibly symmetry protected) 
topologically ordered phases, that are characterized by 
some entangled fixed-point ground state |$f. p .) instead of 
|prod). In this case, the finite range MERA is used to 
coarse-grain the ground state I^gs) into |$f. P ), which is 
then represented by some other tensor network (e.g. an 
MPS, PEPS or scale invariant MERA). Fig. [H;ii) shows 
an example in D=l dimensions. 

Strictly speaking, the exponential decay of correlations 
in the finite range MERA does not accept a simple geo- 
metric interpretation. On the one hand, the ansatz is not 
homogenerous (that is, different tensors describe differ- 
ent length scales) and therefore the relation between the 
distance Dh \(xi,X2) and a correlation function C(xi,X2) 
is no longer be given by Eq. 1131 On the other hand, in 
a finite range MERA all correlations C(x\,X2) for dis- 
tances \xi — X2\ larger than O(2 Z0 ) are zero, as a result of 
the isometric character of tensors in the MERA. This last 
limitation disappears if the tensors at the top layer of the 
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MERA are not constrained to be isometric. 
1 In the MERA for a ground state near a quantum critical 
point, most of the entanglement entropy in a block of size 
L S> £ comes from the bond indices corresponding from 
the first zo = log 2 (£) layers of the finite range MERA, 
which make approximately equal contributions, so that 
the entropy saturates to a constant that depends on the 
correlation length £ as 5'mera(^4) ~ log 2 (£). 



The PEPS resulting from rewriting a scale invariant 
MERA in D > 1 is significantly inhomogeneous (it re- 
quires infinitely many different tensors, one for each value 
of the scale parameter z) and its bond indices have to 
communicate correlations at all length scales. As a result, 
its practical use as a variational ansatz is not clear. 



